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Abstract. Recently, the quaternionic quantum walk was formulated by the first author as a gen¬ 
eralization of discrete-time quantum walks. We treat the right eigenvalue problem of quaternionic 
matrices to analysis the spectra of its transition matrix. The way to obtain all the right eigenval¬ 
ues of a quaternionic matrix is given. From the unitary condition on the transition matrix of the 
quaternionic quantum walk, we deduce some properties about it. Our main results, Theorem 15.31 
determine all the right eigenvalues of a quaternionic quantum walk by use of those of the corre¬ 
sponding weighted matrix. In addition, we give some examples of quaternionic quantum walks and 
their right eigenvalues. 


1 Introduction 

As a quantum version of the classical random walk, the quantum walk has recently been 
studied in various helds. The discrete-time quantum walk in one dimension lattice was 
intensively studied by Ambainis et al. [2] . One of the most striking properties is the spreading 
property of the walker. The standard deviation of the position grows linearly in time, 
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quadratically faster than classical random walk. On the other hand, a walker stays at the 
starting position, i.e., localization occurs. Konno [12] extended the quantum walk to a 
quaternionic quatum walk determined by a unitary matrix whose component is quaternion, 
and presented some properties of the walk. The reviews and books on quantum walks are 
Kempe [I5], Konno im. Venegas-Andraca [30], Cantero et ah [6], Manouchehri and Wang 
ra. Portugal [23] for examples. 

For a general graph, a typical case of discrete-time quantum walks is the Grover walk 
on a graph. The Grover walk on a graph was formulated in m- Emms et ah [H] treated 
spectra of the transition matrix and its positive support of the Grover walk on a graph, and 
showed that the third power of the transition matrix outperforms the graph spectra methods 
in distinguishing strongly regular graphs. Godsil and Guo [10] gave new proofs of the results 
of Emms et ah [S]. 

The discrete-time quantum walk on a graph is closely related to the Ihara zeta function of 
a graph. Zeta functions of graphs originally started from the Ihara zeta function for a regular 
graph by Ihara na. The Ihara zeta function of a graph was studied by many researchers 
H, H, izg, PI, na. Bl. PI. 0. PI- Bass [1] generalized Ihara’s result on the zeta 
function of a regular graph to an irregular graph, and showed that its reciprocal is again 
a polynomial. Inspired by these works, Sato [2S] dehned a new zeta function (the second 
weighted zeta function) of a graph by using not an inhnite product but a determinant. This 
new zeta function and its determinant formula played an essential role in the determination 
of eigenvalues of quantum walk on a graph in Konno and Sato [19]. Ren et al. m found 
an interesting relation between the Ihara zeta function and the Grover walk on a graph, and 
showed that the positive support of the transition matrix of the Grover walk is equal to the 
Perron-Frobenius operator (the edge matrix) related to the Ihara zeta function. Konno and 
Sato [19] gave the characteristic polynomials of the transition matrix of the Grover walk and 
its positive support by using the second weighted zeta function and the Ihara zeta function, 
and so obtained the other proofs of the results on spectra for them by Emms et al. [8]. 

Quaternion, which can be considered as an extension of complex number, was discovered 
by Hamilton in 1843. The eigenvalue problem for matrices whose entries are quaternions 
has been investigated for nearly a century by a number of researchers. A detailed survey 
on the quaternionic matrices can be found in Zhang [3T]. The notable difference from the 
eigenvalue problem for complex matrices is that it is necessary to treat left eigenvalues and 
right eigenvalues separately. Right eigenvalues are well studied by, for example, Brenner [2] 
and Lee [21]. On the contrary, left eigenvalues are less known and not easy to handle as 
commented in Huang and So [I3]. Recently, an extension of quantum walk to the case of 
quaternions was established by Konno [H] in which a model of quaternionic quantum walk 
was suggested and its some properties were presented. In the present paper, we extend a 
discrete-time quantum walk (Grover walk) on a graph to a walk given by a unitary matrix 
whose components are quaternions (quaternionic Grover walk). We derive the condition 
for components of transition matrix of quaternionic Grover walk. Furthermore, we deal 
with right eigenvalue problem, and derive all the right eigenvalues of quaternionic transition 
matrix from those of the corresponding weighted matrix. Our results can be regarded as a 
generalization of [8], [19]. 

The rest of the paper is organized as follows. Section 2 treats the right eigenvalue problem 
of quaternionic matrices. The way to obtain all the right eigenvalues of a quaternionic matrix 
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is given. In Section 3, we provide a summary of the Ihara zeta function and the second 
weighted zeta function of a graph, and present their determinant expressions. In Section 
4, we show a brief account of the discrete-time quantum walk on a graph, and dehne a 
generalization of it to the case of quaternions, which we call the quaternionic quantum walk 
on a graph. From the unitary condition on the transition matrix, we deduce some properties 
about it. In Section 5, our main results. Theorem 15.31 determine all the right eigenvalues 
of a quaternionic quantum walk by use of those of the corresponding weighted matrix. In 
addition, we give some examples of quaternionic quantum walks and their right eigenvalues. 

2 Right eigenvalues of a quaternionic matrix 

In this section, we shall show the way to calculate all the right eigenvalues of a quaternionic 
matrix. Since quaternions do not mutually commute in general, we must treat left eigenvalue 
and right eigenvalue separately. Since M”v = vA” (n > 1) holds for an right eigenvalue A of 
M and an eigenvector v corresponding to A, right eigenvalue is more suitable to calculate the 
stationary measure of the quaternionic quantum walk than left eigenvalue. Therefore, we 
treat only right eigenvalues in this paper. Although we cannot apply the ordinary character¬ 
istic polynomial for a quaternionic matrix, we can calculate all the complex right eigenvalues 
of it by embedding it in a complex matrix twice the size of it. From these eigenvalues, we can 
obtain all the right eigenvalues of the quaternionic matrix. An exposition of these contents 
can be found in 0- Ell gives an overview of quaternionic matrix theory. 

Let El be the set of quaternions. El is a noncommutative associative algebra over M, 
whose underlying real vector space has dimension 4 with a basis which satisfy the 

following relations: 

= k‘^ = —1, ij = —ji = k, jk = —kj = i, ki = —ik = j. 

For X = Xq + Xii + X 2 j + x^k G H, we call Xq the real part, Xq + Xii the complex part, xii + 
X 2 j + X 3 k the imaginary part oi X, and set Re x = Xq, Co x = Xg + x^i, Im x = Xii + X 2 j + X 3 k 
respectively, x* denotes the conjugate of x in El which is dehned by x* = xq — xii — X 2 j — xgk. 
We call |x| = \/xx* = y/x*x = \/xg + xf -|- x| -|- x| the norm of x. Indeed, | • | satishes 

(1) |x| > 0, and moreover |x| = 0 <;=> x = 0, 

(2) \xy\ = \x\\y\, 

(3) \x + y\< |x| + \y\. 

For a nonzero element x G El, x~^ = x*/|xp. Hence, El constitutes a skew held. 

We can present arbitrary quaternion x by two complex numbers x = a+jh uniquely. Such 
a presentation is called symplectic decomposition. Two complex numbers a and b are called 
simplex part and perplex part of x respectively. Symplectic decomposition is also valid for 
a quaternionic matrix, namely a matrix whose entries are quaternions. Mat(mx?7,,EI) de¬ 
notes the set of mxn quaternionic matrices and Mat(n,H) the set of nxn quaternionic 
square matrices. For M G Mat (mxn, El), we can write M = A -|- jB uniquely where 
A, B G Mat (mxn, C). A and B are called simplex part and perplex part of M respectively. 
We dehne 'ijj to be the map from Mat (mxn, El) to Mat(2mx2n, C) as follows: 

A —B^ 

: Mat(mxn,H)—>■ Mat(2mx2n,C) M i—)■ , 
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where A is the complex conjugate of A. Then tjj is an M-linear map. We also have: 
Lemma 2.1 Let M G Mat(mxn, H) andN G Mat(nxm, H). Then 


V'(MN) = ?/^(M)V'(N). 

Proof. Let M = A + jB and N = C + jD be symplectic decompositions of M and N. 
Then, 

MN = (A + jB)(C + jD) = AC + AjD + jBC + jBjB. 

Since Xj = jX for arbitrary complex matrix X, it follows that: 

MN = AC - BD + j(AD + BC), 


and therefore 


On the other hand, 


V;(MN) 


AC-BD -AD-BC 
AD + BC AC-BD 




A 

-B 


C 

-D 


AC - BD 

-AD - BC 

B 

A 


D 

C 


BC +AD 

-BD + AC 


Thus = ?/;(M)V’(N) holds. □ 


Proposition 2.2 If m = n, then i/j is an injective W-algebra homomorphism. 

Proof. By Lemma [2.11 is an M-algebra homomorphism. Injectivity of if is clear. □ 
For a quaternionic square matrix M = (M^) G Mat(n, H), the conjugate M* = ((M*)jj) 
is dehned by We notice that '0(M*) = '0(M)*, where the right-hand side 

denotes the conjugate transpose of the complex matrix. A quaternionic square matrix M is 
said to be unitary if M*M = MM* = I. Since if is an injective M-algebra homomorphism, 
M*M = I implies MM* = /. 

Remark 2.3 '0(^lM) = ^^(M) does not hold in general, where ^TM denotes the transpose of 

M. 


Let El”' be a right vector space, and M G Mat(n, H), v G H”, A G El satisfy Mv = vA. A 
is said to be a right eigenvalue of M, and v a right eigenvector corresponding to A. We denote 
by crr(M) the set of all the right eigenvalues of M. We have M(vg) = vAg = vg(g“^Ag) 
for any g G H* = El — {0}. Hence, the conjugate class A®'* = {g“^Ag | g G H*} is contained 
in (Tr(M). We notice that if A G M, then A®** = {A}. Since |g“^xg| = |x|, the group 
homomorphism: 

p : H* —)■ Ant El, p{q)x = qxq~^ 

gives orthogonal transformations p(g) (g G H*) on El. Let Hr = M and Hp = Mi©Mj©Mfc. 
Then H = Hp©Hp and Hr T Hr. We call the elements of Hr pure quaternion. One can see 
that p{q)x = X {x E Hr) and p(g)HR = Hr, hence p also gives orthogonal transformations 
p(g) (g G H*) on Hr. Moreover it is well-known that det(p(g)) = 1 and that the group 
homomorphism: 

pIhp : H* ^ 50(Hr) = 50(3) 
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is surjective (for details, see for example [I],[7]). Therefore, H* gives all the rotations in 

Hp = MT 

If A G El — R, then we may write A = Aq + A' (Aq G Hp, A' G Hp) with A' 7 ^ 0. Then 
p{q)X = qXq~^ = Aq + qX'q~^. It follows that A'®* is a 2-sphere centering at the origin in Hp 
since the following equations hold: 


A®* = { q -^ X'q I g G H*} = p(H*)A' = 50(Hp)A'. (2.1) 

Thus the intersection of A®* and the complex axis consists of just two points, A®*nRi = 
{±|A'| 2 }. To summarize our discussion, we hnd out the following: 

Proposition 2.4 Let M g Mat(?7,,H). The M*-conjugate class of any right eigenvalue of 
M which is not real is contained in (Jr(M), and has exactly two complex numbers that are 
conjugate. 

We shall show that one can calculate all the complex right eigenvalues of a quaternionic 
matrix by use of ordinary determinant. 


Proposition 2.5 Let M g Mat(?7,,H) andp G C. Then 


Mv = \p ^ '^(M) 


u 


u 

w 


w 


p where v = u-|-jw (u, w G C"'). 


Proof. Let M = A -|- jB be the symplectic decomposition. Then we have: 

Mv = (A -|- jB) (u -|- jw) = (Au — Bw) -|- j (Bu -|- Aw). 


Hence Mv = vp is equivalent to the following equation: 


A 

-B 


u 


u 

B 

A 


w 


w 


□ 

In the same way as the ordinary triangularization of a complex matrix, the triangular- 
ization of a quaternionic matrix can also be obtained. Moreover, we can triangularize a 
quaternionic matrix into a triangular matrix whose diagonals are complex as follows. The 
detail of the proof can be found in I2I1- 


Lemma 2.6 For arbitrary M. G Mat(n, H), there exist V G Mat(?7,,H) withY*Y = VV* = 
In and Ai,..., A„ G C such that 


V*MV 


Ai * 

0 A 2 


0 0 ■ ■ ■ Xn 


( 2 . 2 ) 


As a consequence of Proposition 12.51 and Lemma 12.61 we can deduce the following: 
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Proposition 2.7 Let M g Mat(n,H). Then 2n eigenvalues o/V’(M) G Mat(2n,C) occur 
in conjugate pairs: 

•^ 1 , -^ 1 ; -^ 2 ; -^ 2 ; • • •, ^ni ^ni 

counted with multiplicity. 


Proof. The proof can be found in m- Let T be the right hand side of 02.21) and 
T = T'^ + its symplectic decomposition. Then, 

det(Al2„ - V'(M)) = det(Al2n - V'(V)-V(M)V'(V)) 

= det(Al2n - ^{'V*MV)) 


* 

* 


= det(Al 2 n - '0(T)) = det 


Alon - T‘ 


A-Ai 
0 


* 

A — A 2 


* 

* 


rpp 

Al 2 n-T^. 

0 * 

0 0 


□ 


= det 


0 

0 

0 


0 

* 

0 


A — A^ 0 0 

* A — Ai * 

* 0 A — A 2 


0 

* 

* 


0 0 0 0 
= (A-Ai)(A-AT)---(A-A„)(A-^). 


0 ■ • • X — Xn 


Hence, diagonal elements in T and their complex conjugates are eigenvalues of ■^(M). 
We notice that we may replace Aj by its complex conjugate in 02 .21) because of Proposition 
12.41 Combining Proposition 12.4112.5112.71 we can state about eigenvalues of a quaternionic 
matrix as follows: 


Theorem 2.8 For arbitrary guaternionic matrix M G Mat(n, H), 2n complex right eigen¬ 
values Ai,..., A„, Ai, ..Xn o/M exist counted with multiplicity, which can he calculated by 
solving det(Al 2 n — V’(M)) = 0. The set of right eigenvalues cTr(M) is given by (Tr(M) = 
Afu---UAf. 


Example 2.9 M 


0 1 

-1 0 


det(Al 4 — V’(M)) = det 

Example 2.10 M = 

det(Al 4 — -^(M)) = det 


A -1 
1 A 
0 0 
_0 0 

1 O' 

0 i 
' "A - 1 
0 
0 
0 


0 0 
0 0 
A -1 
1 A 


(A-i)2(A + i)2. a^(M)=iH‘. 


0 0 0 

X-i 0 0 

0 A-1 0 

0 0 A + i 


(A-1)2(A-*)(A + *). 
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a,(M) = 


Example 2.11 M= ^ ^ ^ ^ 

j OJ [0 OJ 

[A -i 0 O' 

det(Al 4 — '^(M)) = del n \ ^ = 

0 0 A « 

-1 0 0 A 








(A2-*)(A2 + *) 



3 The Ihara zeta function of a graph 

Let G = iy (G), -E(G)) be a finite connected graph with the set V (G) of vertices and the set 
E{G) of unoriented edges uv joining two vertices u and v. We assume that G has neither 
loops nor multiple edges throughout. For uv G E{G), an arc {u, v) is the oriented edge from u 
to V. Let D{G) = { {u,v), {v,u) \ uv G E{G)} and |E(G)| = n, |E(G)| = m, |iA(G)| = 2m. 
For e = (m, v) G E(G), o(e) = u denotes the origin and f(e) = v the terminaloi e respectively. 
Furthermore, let e~^ = {v,u) be the inverse of e = {u,v). The degree degn = degc v of a 
vertex n of G is the number of edges incident to v. A path P of length £ in G is a sequence 
P = {ei, ■■■ ,ei) of i arcs such that Cj G D{G) and t^ei) = o(ei+i) for i G {1, •••,£ — 1}. 
We set o(P) = o(ei) and t{P) = t^ef). |P| denotes the length of P. We say that a 
path P = (ci, • • • , e^) has a backtracking if e^+i = e~^ for some i{l<i<l — 1), and that 
P = (ei, • • • , e^) has a tail if ei = ef^. A path P is said to be a cycle if t{P) = o(P). The 
inverse of a path P = (ei, • • • , e^) is the path {ef^, • • • , ef^) and is denoted by P~^. 

We introduce an equivalence relation between cycles. Two cycles Gi = (ci, • • • , e^) and 
C 2 = (/i, • ■ ■ ) Z^) are said to be equivalent if there exists k such that fj = for all j 
where indices are treated modulo Let [G] be the equivalence class which contains the 
cycle G. Let P*” be the cycle obtained by going r times around a cycle B. Such a cycle is 
called a power of B. A cycle G is said to be reduced if G has no backtracking and no tail. 
Furthermore, a cycle G is said to be prime if it is not a power of a strictly smaller cycle. 

The Ihara zeta function of a graph G is a function of t G C with \t\ sufficiently small, 
defined by 

c 

where [G] runs over all equivalence classes of prime, reduced cycles of G. 

Let G be a connected graph with n vertices and m unoriented edges. Two 2m x 2m 
matrices B = (Be/)e,/ez)(G) and Jq = (Je/)e,/eD(G) ai’e defined as follows: 


Be/ - 


1 if t(e) = o(/), J 
0 otherwise, 


1 iff = e-\ 
0 otherwise. 


Then the matrix B — Jq is called the edge matrix of G. 
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Theorem 3.1 (Hashimoto ||I2]; Bass [i]) Let G be a connected graph. Then the recip¬ 
rocal of the Ihara zeta function of G is given by 

Z(G, t)-' = det(l 2 ^ - t(B - Jo)) = (1 - t^y-^ det(I„ - tA + - !„)), 

where r and A are the Betti number and the adjacency matrix of G, respectively, and D = 
(Du^)„^„gy(G) is the diagonal matrix with = degw for arbitrary u G V{G). 

Consider an n x n complex matrix W = {yVuv)u,v£V{G) with (n, n)-entry eqnals 0 if 
{u,v) ^ L){G). We call W an weighted matrix of G. Fnrthermore, let w{u,v) = for 
u,v E V (G) and w{e) = w{u, n) if e = (m, v) G D{G). 

For an weighted matrix W of G, a 2m x 2m complex matrix B^ = is 

defined as follows: 

gW ^ f w{f) if t{e) = o{f), 

[0 otherwise. 

Then the second weighted zeta function of G is defined by 

Zi{G,w,t) = det(l 2 m - - Jo))"^ 

If w{e) = 1 for any e G D{G), then the second weighted zeta fnnction of G coincides with 
the Ihara zeta fnnction of G. 

Theorem 3.2 (Sato (25j ) Let G be a connected graph, and let W be a weighted matrix of 
G. Then the reciprocal of the second weighted zeta function of G is given by 

ZyG,w,t)-^ = (l-f2)™-"det(I„-fW + f2(D^-!„)), 

where n = |B(G)|, m = |F'(G')| and D^, = is the diagonal matrix with = 

w{e) for arbitrary u G V{G). 

e\o{e)=u 

Remark 3.3 For later use, we mention that taking transpose, the following eguation also 
holds: 

det(l 2 ^ - f(^B^ - Jo)) = (1 - det(I„ - fW + - R)). 

4 The transition matrix of a quaternionic quantum 
walk on a graph 

A discrete-time quantnm walk is a qnantnm process on a graph whose state vector is governed 
by a nnitary matrix called the transition matrix. An important example of the qnantnm 
walk on a graph is the Grover walk, which was introdnced in m- Let G be a connected 
graph with n vertices and m edges. Set du = degn for u G V{G). The transition matrix of 
the Grover walk = {\]^f°)ej&D{G) of G is defined by 

f 2/(i,(/)(= 2/d„(e)) if«(/) = o(e) and / 7 ^ 

Uf/“=1 2/d„/)-l it/ = e-‘, 

I 0 otherwise. 


U is called the Grover matrix. We denote by Spec{A) the multiset of spectra of a complex 
square matrix A counted with multiplicity. Let T = {T^uv)u,v^v{g) be the n x n matrix 
dehned as follows: 

^ f 1/du if {u,v) e D{G), 

0 otherwise. 

In |8], Emms et ah determined the spectra of by use of those of T. 

Theorem 4.1 (Emms, Hancock, Severini and Wilson [S]) LetG be a connected graph 
with n vertices and m edges. The transition matrix has 2n eigenvalues of the form: 

A = At i 

where At is an eigenvalue of the matrix T. The remaining 2{m — n) eigenvalues are 

±1 with egual multiplicities. 

Now, we extend the Grover walk to the case of quaternions. A discrete-time quater- 
nionic quantum walk is a quantum process on a graph whose state vector, whose entries are 
quaternions, is governed by a quaternionic unitary matrix called the quaternionic transition 
matrix. Let G be a hnite connected graph with n vertices and m edges. The guaternionic 
transition matrix U = (Ue/)e,/GD(G) of G is dehned by 

f <?(e) if t{f) = o(e) and / ^ e"^ 

Ue/ = < q{e) - 1 if / = e“^ 

y 0 otherwise, 

where g is a map from D{G) to H. The unitary condition on U is equivalent to the following 
equations: 

Lemma 4.2 Let g(e) = qo{e) -I- qi{e)i + q 2 {e)j + qz{e)k. Then U is unitary if and only if the 
following five equations hold: 

qo{eY + qi{ef + q 2 {,ef + q^i^ef - = 0, (4.3) 

do{e) 

(?o(e)go(/) + qi(,e)qi{f) + g 2 (e)g 2 (/) + g3(e)g3(/))do(e) - (go(e) + go(/)) = 0, (4.4) 

{-qo{e)qi{f) + qi{e)qo{f) - q 2 {e)qz{f) + q3{e)q2{f))do{e) - {qi{e) - qi{f)) = 0, (4.5) 
(-go(e)g2(/) + g2(e)go(/) - g3(e)gi(/) + gi(e)g3(/))do(e) - (g2(e) - g2(/)) = 0, (4.6) 
{-qo{e)q3{f) + q3{e)qo{f) - qi{e)q2{f) + g2(e)gi(/))do(e) - {q3{e) - q3{f)) = 0, (4.7) 

where e and f are different arcs and have the same origin. 

Proof. Observing that 

(UU-)./ = 

g&D{G),t{g)=o(e)=o{f) 

one can readily check by direct calculations that (UU*)ee = 1 is equivalent to 04.31) . and 
(UU*)e/ = 0 to g3D,(!ll5]),(IMD and gT]). If o(e) ^ o(/), then clearly (UU*)e/ = 0 holds. 

□ 
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From fl4.3p . it follows that qo{e) must satisfy 


0 < go(e) < 


2 


d‘o{e) 


Theorem 4.3 If U is unitary, then q{e) = q{f) holds for any two arcs e, / G D{G) with 
o(e) = o{f). 

Proof. We set qo{e) = x/du, qo{f) = y/du{0 <x,y < 2) and q(e) = ^{qi{e), q 2 {e), qsie)), 
where u = o{e) = o{f). Then fl4.4p turns into: 

(x-l)fa-l)-l _ 0 (4 8) 

On the other hand, fl4.5l) . fl4.6l) . fl4.7p can be expressed all at once as follows: 

{y - l)q(e) -{x- l)q(/) - d„q(e)xq(/) = 0. (4.9) 

We show the assertion divided into three cases depending upon q(e) and q(/). 

If q(e) = q(/) = 0, then fl4.9p trivially holds and fl4.8p turns into (a: — 1 )(|/ — 1) — 1 = 0. 
In addition, by fl4.3p we have 


Thus, {x,y) must be (0,0) or (2,2). 

If only either one of q(e) or q(/) equals 0, say q(/) = 0, then fl4.9p becomes (?/ —l)q(e) = 
0. Hence y = I and q{f) = l/du- But from 04.81) it follows that —l/du = 0, thereby reaching 
a contradiction. 

If both q(e) and q(/) do not equal 0, then by 04.91) it follows that: 


{y - l)q(e) -{x- l)q(/) = 0, (4.10) 

q(e)xq(/) = 0. (4.11) 

By 04.111) . q(e) and q(/) are parallel, so we may write q(/) = cq(e) for some nonzero real 
number c. Then O4.10p implies y — 1 = c(x — 1), and 04.8p yields: 

c{{x-l)^ + dj\q{e)\^} = l (4.12) 

Finally, substituting x = qo{e)du for 04.121) and using 04.31) . we get c = 1. □ 

Conversely, if g(e) = g(/) for any two arcs e, / G D{G) with o(e) = o(/), then 04.4p turns 
into 04.31) . and moreover, 04.51) . 04.Op . 04.7p hold trivially. Thus we obtain: 


Corollary 4.4 U is unitary if and only if 

qoief + gi(e)2 + q2{ef + qs^ef - = 0, 

«o(e) 

and g(e) = g(/) for any two arcs e, / G D{G) with o(e) = o(/). 
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Hereafter we impose the following condition on q{e): 

q{e) does not depend on u. (4-13) 

e:o(e)=u 

Hence, we may write a = ao + aii + a 2 j + ask = J2e-o{e)=u^(^) arbitrary u G V{G). 
Theorem 14.31 and fl4.13p imply g(e) = a/do(e)- If fhe imaginary part of a equals 0, then U 
must equal either —Jq or Therefore, U can be regarded as a quaternionic extension 

of the Grover walk. 


5 Eigenvalues of the transition matrix of a quater¬ 
nionic quantum walks on a graph 


We shall calculate all the right eigenvalues of U. By Theorem 12.81 it suffices to solve 
det(Al 4 m — "^(U)) = 0. Consider the set of quaternionic conjugates 

U^‘ = {q-^Vq I q G H*}, 

of U. By fl2.ip and the subsequent argument, a®*nC = a± where a± = ao±\/a? + 0 ^ + a\ i. 
Let U± = {(U±)ef)e,feD{G) be two complex matrices defined as follows: 


(U±)e/ = < 


a± 

do{e) 

a± 

do{e) 

0 


- 1 


if f(/) = o(e) and / 7 ^ e ^ 

if / = e-^ 

otherwise. 


We notice that U_ = U+. Then, U'®*n Mat(2m, C) = {U+, U_}. Therefore, for the 
quaternion q such that q~^aq = «+, we have: 

det(Al4m - ^/'(U)) = det('0(gl2m)"^(Al4m - ^/J{\J))'^/J{ql 2 m)) 

= det(Al4m - '0(gl2m)"V(U)'0(gl2m)) 

= det(Al4m - 
= det(Al4m - i>{\J+)). 


Thus, under the condition fl4.13p . we can calculate all the right eigenvalues of U by calculating 
eigenvalues of V’(U+). However, 


det(Al4m - V’(U+)) 


det ^ 

Al2m O 

O Al2m 

— 


det(Al 2 m — U+) det(AI 

2m -U_). 


(5.14) 


Hence, S'pec('0(U)) = S'pec('0(U+)) = S'pec(U+)US'pec(U_). Since S'pec(U_) is obtained 
from Spec{U+) by conjugating their elements, S'pec('0(U)) consists of eigenvalues of U+ and 
their complex conjugates. Thus by Theorem 12.81 we obtain: 
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Theorem 5.1 Let \J is a quaternionic transition matrix on a graph which satisfies ( 4-i3 ). 
Then crr(U) is given by: 

a,(U) = Afu---UAf, 
where Ai, • • •, A„ are eigenvalues of U+. 

Now we apply Theorem 13.21 to the eigenvalue problem for U. Let us dehne the nxn 
complex matrix Wj. = as follows: 

av 1 J ^ if («,!.) 6 C(G). 

0 otherwise. 


Proposition 5.2 Let G be a connected graph with n vertices and m edges. Then, for the 
quaternionic transition matrix U of G, 

detiXhm - U±) = (A^ - l)”^-"det((A2 + a± - 1)1, - A^W±). 


Proof. Let G be a connected graph with n vertices and m edges. Corresponding to W±, 
we dehne B„± = {(B±'^)ef)e,f&D{G) ^ Mat (2m, C) to be as follows: 




a± 


do{f) 

0 


if f(e) = o(/), 
otherwise. 


We readily see that U± = — Jo- Then, by Theorem 13.21 and the subsequent remark, it 

follows that: 

det(l2m - tlJ±) = det(l2r„ - - Jq)) 

= (1 _ ^ 

= (1 _ t^)rn-n ^ 

= (1 - det((l + (a± - l)t^)In - t^W±). 

Now, let f = 1/A. Then we have 

det (l 2 m - “ i) ((^ + ~ ■ 

This implies the following conclusion: 

det(AI2,^ - U±) = (A^ - 1)"^-" det((A2 + a± - 1)1, - A^W±). 


□ 

Using Proposition 15.21 we can calculate the spectra of "(/(U) as follows. 

Theorem 5.3 Let G be a connected graph with n vertices and m edges. S'pec('0(U)) has 
4m elements. If G is not a tree, then 4n eigenvalues of them are of the form: 

fi± \/fjf — 4(a+ — 1) V ± \/u‘^ — 4(a_ — 1) 
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where /i G S'pec(), i/ G Spec(^"W_). The remaining 4(m — n) eigenvalues o/'0(U) are 
±1 with egual multiplicities. If G is a tree, then S'pec('^(U)) is as follows: 


Specitl^fU)) 

^ ^jj? — 4(a+ — 1) 1 / ± -^1/2 — 4(a_ — 1) 

^ I 2 ’ 2 

-{ 1 , 1 ,- 1 ,- 1 }. 


/iG^pec(^W+), i/G^pecfW_))| 


Proof. By fl5.14p and Proposition 15.21 we have 


det(AI4„^-^(U)) = (A2-1)2™-2- Yl {X^ + a+-l-pX) JJ (A" + a_-l-i/A). 

/i£5pec(^^W+) u£Spec(^"W-) 

Solving A^ + a+ — 1 — /rA = 0 and A^ + O- — 1 — z/A = 0, we obtain 

^p? — 4(a+ — 1) n ± — 4(a_ — 1) 

~ 2 ’ 2 ' 

If G is not a tree, then m >n and hence the statement clearly holds. If G is a tree, then 
m = n — 1. Since det(Al4m — V’(U)) is a polynomial of A, npe 5 pec(rw+)(-^^ + a+ — 1 — 
tX) niye 5 pec(Tw_)(-^^ + a_ — 1 — z/A) must have the factor (A^ — 1)^. Thus the statement 
holds for trees. □ 


Remark 5.4 Since Specff^ J) is obtained from S'pec(^W+) by conjugating their elements, 
it suffices to calculate only the eigenvalues of Specif'X^^). 


Example 5.5 Let G = K^, the complete graph with 3 vertices, and a = l + l/2f + \/2/2j — 
l/2k. We give a weighted matrix as follows: 


W = 


a 

2 

a 

L2 


a 

2 

0 

a 

2 


Then the guaternionic transition matrix U is as follows: 


U = - Jn = 


0 

- - 1 

a 

2 

0 

0 

0 


a 

2 


0 

0 

0 

0 


0 

0 

0 

- -1 
Q 
2 
0 


- -1 
0 
0 
0 


a 

2 

0 

0 

0 

0 

- -1 
2 


0 

0 

0 


2 

- -1 

2 


0 


Then a± = l±f and 


Spec{^W+) = 




, . 11.11 
- + -i,-- + -I 


13 


















and 


Spe.(U.) ^ {1.. j 

Spec(U_) ^ + + 

Thus we obtain: 

,.,> ,, H. (-i + vi -i-viT (-1-V7 -i + viT 

MU) = {i}u"u(^^^ + ^^.j + . 

Next, we state speetra of the Grover matrix. Let a = 2. Then U is the Grover matrix 


0 0 0 0 1 0 

0 0 0 1 0 0 

Gro^ 1 0 0 0 0 0 

0 0 0 0 0 1 ’ 

0 0 1 0 0 0 

_0 1 0 0 0 0 _ 

Furthermore, the transition matrix T of the simple random walk of is given as follows: 

To i il 

T = 1 0 i 

i i 0 

L2 2 


Then 


Spec{T) = ll,~~ 


and S'pec(U*^™) = <1,1, 


-l±^/3i -l±V3i 


Thus we obtain: 


ariV^n = {1}U 


-1 + V3i 


Example 5.6 Let 0 = 84 , the star graph with 4 vertiees, and a = 4/3+l/3i+2/3j+\/3/3A;. 
We give a weighted matrix as follows: 


W = 


0 0 0 a 

0 0 0 a 

0 0 0 a 

a a a n 

.3 3 3 


Then the guaternionie transition matrix U is as follows: 


- - 1 

3 


U = 'B,„ - Jn = 


a — 1 
0 
0 
0 
0 
0 


a _ -| 

3 ^ 

0 


0 

0 

a — 1 
0 
0 
0 


0 

0 

0 

0 

a — 1 
0 


14 



Then a± = 4/3±2\/2/3i and 


r. X ) 4 2^2 4 2^2 1 ^ X f 4 2^2 4 2^2 

^pec(^W+) = 0,0,- + ^—%,-- - —I L Spec(^W.) = < 0,0, 


3 3 3 3 


3 3 *’“3 + ^' 


and 

Spec{V+) = <j ±( 


1 V2. 


1 y ^\ 1 2 ^ 2 . 1 2^/2 . 


\ / i A i 

x/3 x/3V’=^i7!“7!7'^’3 + 


3 3 


^,-u-{l,-l}. 


O /TT ^ ; , / 1 , M . 1 2^2. 1 , 2^2. ^ ^ 

Spec(U_) = ,( ±( — + _,j, ±(^_ + _,j, 1, - - -1 ^ - {1, -1}, 


T/itis we obtain: 


nv: I 1 ( 1 ^X' A 272.', , 

= 17!■ 7!V A"7! + 7!V A5 + —4 A 3 3 


1 2^2 


Next, we state speetra of the Grover matrix. Let a = 2 and U = Then is 


given as follows: 


UGro ^ 


0 
-1 
3 

0 0 


2 

L 3 


1 

0 


0 0 0 0 


0 


2 

3 

0 

-1 

3 


0 


0 

1 

0 


0 0 0 0 


0 


2 

3 

0 

-1 

3 


^ 0 
3 ^ 

0 0 


Furthermore, the transition matrix T of the simple random walk of S 4 is 


Then 


0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

1 

.3 

1 

3 

1 

3 

0 


Spec{T) 


0 , 0 ,- 1,1 


and S'pec(U^™) 


ih i, ±i, 1, —1 


Thus we obtain: 

= {1}U{-1}UF. 
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